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Abstract 

O | We analyze the properties of asymptotically AdS electrically charged black brane solutions in a 

consistent truncation of the Af = 4 + , D — 5 Romans' gauged supergravity which contains gravity, SU(2) 
*Jm and U(l) gauge fields, and a dilaton possessing a nontrivial potential approaching a constant negative 

value at infinity. We find that the Z7(l) X U(l) solutions become unstable to forming non-Abelian hair. 
These configurations emerge as zero modes of the Abelian solutions at critical temperature and a critical 
(nonvanishing) ratio of the electric charges and can be viewed as holographic p— wave superfluids. 



1 Introduction 



Recently, considerable effort has been put into extending AdS/CFT correspondence beyond high-energy 
physics by constructing gravity models that are conjectured to be dual to various condensed matter systems. 
This has lead to the discovery of holographic superconductors and holographic superfluids, describing con- 
densed phases of strongly coupled, planar, gauge theories. Studying such models involves the construction 
of electrically charged black holes in an asymptotically AdS spacetime, which, below a critical temperature 
become unstable to forming hair. That is, a phase transition occurs to a superconductor/superfluid state, 
in which a sufficiently large U(l) charge density triggers the spontaneously breaking of the U(l) symmetry. 
Then an operator charged under the U(l) acquires a nonzero expectation value (see e.g. [1] for a review of 
these aspects). 

For p— wave superconducting black holes, the condensing operator is a vector and hence rotational sym- 
metry is broken. Such black hole solutions have been constructed using either charged non-Abelian vector 
fields [2] or, alternatively, charged two- forms [5]. However, most of the studies in the literature have assumed 
an ad hoc construction of the lagrangean of the gravitational system, without a clear connection with a given 
supergravity model, which makes it rather difficult to describe precisely the application of the AdS/CFT 
dictionary. 

At the same time, the gauged supergravity models generically contain non-Abelian vector fields, which 
may suggest the existence of p— wave superconducting black hole solutions. The case of Af = 8, D = 5 
gauged supergravity [3], [3] is of particular interest, given its connection with Af = 4 U(N) super- Yang-Mills 
theory in 3 + 1 dimensions. The bosonic sector of this theory consists of the metric, twenty scalars and fifteen 
50(6) Yang-Mills (YM) gauge fields^. Solutions of Af = 8, D = 5 model have been considered by several 



1 Note that the field content of the full Af = 8, D = 5 gauged supergravity is richer. However, a number of bosonic fields 
can be consistently set to zero [5]. 
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authors for various consistent truncations, with subgroups of SO(6) (see e.g. [B] and the references therein). 
However, to our knowledge, to date no attempt has been made to construct non-Abelian superconducting 
black hole solutions in this context. 

This paper is aimed as a first step in this direction, by taking a consistent truncation of the Af = 8 
model corresponding to Af = 4 + , SU(2) x U(l) Romans' gauged supergravity, with a single scalar field <j> 
possessing a potential V(cj>) which is the sum of two Liouville terms. The scalar (f> approaches asymptotically 
a constant value 4>q corresponding to an extremum of the potential, dV/d<f>\^ — 0, which yields an effective 

cosmological constant A e // = 2V(4>o ) < 0- It turns out that the basic properties of the Af — 4 + solutions 
with non-Abelian fields are rather similar to those found for pure D — 5 Einstein- YM- A system [8], [9]. In 
particular, we find evidence for the existence, at low temperatures, of a superfluid state with a normalisable 
non-Abelian condensate. 

Since Romans' theory arises as a consistent Kaluza-Klein truncation of the type IIB supergravity on 
an S 5 [10] and as a consistent compactification of D — 11 supergravity [TTJ, this shows the existence of 
holographic superfluids in D = 10, 11 supergravities. 



2 The J\f = 4 + , D = 5 Romans' gauged supergravity 

The bosonic sector of the Af = 4, D = 5 Romans' gauged supergravity [7] consists of gravity, a scalar <j>, 
an SU(2) YM potential (with field strength F^J = d^A^ - d u A^ ] + g YM e IJK A { U J) A^ ] and g Y M the 
SU(2) gauge coupling constant), an Abelian potential (f^ v — d^B u — d v B^ being the corresponding field 
strength), and a pair of two-form fields. These two form fields can consistently be set to zero, which yields 
the bosonic part of the action 

hulk = ^J m d 5 x^—g(±R - \d^d^ - l e ^ F $F^ v - je" 4 ^/^/^ (2.1) 

--l= e ^P°T F $F£B T - V(<j>) 



where a = ^J\. Here V(4>) = -\g YM (e^ 2 ^ + 2 v / 2^-e a *J is the dilaton potential, g M being the U(l) 
gauge coupling constant. 

As discussed in [7], this theory has three canonical forms, corresponding to different choices of the gauge 
coupling constant gM- The case of interest here corresponds to the Af = 4 + version, in which gM = gYAi/V^ 
and thus the dilaton potential is 

V{<t)) = ~ l -g YM (e-^ + 2^). (2.2) 

The field equations are obtained by varying the action (|2.1[) with respect to the field variables g^ v , A^P , £> M 
and 4> 

R„u - \guvR = 2 V, V 2 - ^ F (J0 F (D^ + ae'^f^r ~ ^ = 0, (2.3) 

OJe-^ ^e^P^FyjF^l) = 0, DJe^F^") ]=e» Vf " TT F^ f aT = 0, 

where the energy-momentum tensor is defined by 

Tpv = d^d^-^g^Md^-g^V^) (2.4) 
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The scalar potential has exactly one extremum at <f> = 0, corresponding to an the effective cosmological 
constant A. e ff = — w = 2V(0) = —j9ym- Then the effective AdS length scale is fixed by the non-Abelian 
gauge coupling constant, I = 2v2/ 9ym- 

As usual, one supplements (|2.ip with a boundary term 

' 1 d 4 xV^hK-^-[ d 4 xV^h (]w ((/)) + ^-TZ) , (2.5) 



8?r JdM " 8lT JdM " 4 

where apart from the Hawking- Gibbons surface term we include also a counterterm part which is required 
to regularize the total action and the global charges. In the above relation, 1Z is the Ricci scalar for the 
induced metric h of the boundary, K is the trace of the extrinsic curvature, while W(4>) = e 2a ^ + 2e~ a< ^ (this 
expression of the counterterm was derived in [12] , in a more general context). 

Then, as in the well known pure- AdS case [13] , one can construct a divergence-free boundary stress tensor 
Tjj from the total action I=/buik+^bound by defining 

T« = = -^( K a ~ Kh a ~ \**W{<t>) + \ E al ( 2 -6) 

where is the Einstein tensor of the boundary metric, K irj = — 1 /2{V{rij + V^n,) is the extrinsic curvature, 
with ri 1 being an outward pointing normal vector to the boundary. 
Thus, a conserved charge 

5 = £d 3 S a e"T ab , (2.7) 

can be associated with a surface E (with normal n a ), provided the boundary geometry has an isometry 
generated by a Killing vector £ a . For example, if £ = d/dt is a timelike Killing vector, then 0^ is the 
conserved mass M.. 



3 The uncondensed phase 
3.1 The solutions 

We start with a discussion of the basic properties of the Abelian black brane solutions of the Af = 4 + 
Romans' model. They can be found as a particular limit of the black holes obtained in [14] in the so-called 
STU model. In the general case these black holes possess three different U(l) charges and two independent 
scalars. After setting one scalar to zero and taking two gauge fields to be equal, one finds after a suitable 
field redefinition, the following black brane solution of the eqs. 



ds 2 = H(r) 1 / 3 (^-+ r 2 {dx 2 + dy 2 + dz 2 )) - H{r)~^ f{r)dt 2 , (3.8) 



with 

20 2 4o 2 M 1 

H(r)=H 2 (r)K(r), ff( r ) = l+J^, #( r ) = i + -JL-, /( r ) = -_ + M r 2 U(r) (3.!» 



and the matter fields 



(r) = log , B = B t (r)dt, = A t (r)S I3 dt, (3.10) 

with B t (r) = ^ - - 2 M l r 2 , A t (r) = 



Mr 2 ' v ; 2Q 2 + Mr 2 ' 

This solution is written in terms of three parameters (M,Q,q), corresponding (up to some factors) to the 
global mass and two electric charges. 
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In what follows, to avoid cluttering our expressions with complicated factors of qym, we use the observa- 
tion that the above solution is left invariant by the transformation r — ¥ Xr, gyM — ► 9ym /A, (q, Q) — > X(q, Q) 
and (x, y, z) — > X(x, y, z), and we set gyu = 1 without any loss of generality. 

The horizon is located r = ru, with th the largest positive root of the equation /(r) = 0, which reduces 
to r% + 4(2-±^)rff + ^(Q 4 + Aq 2 Q 2 - 2M 3 )r% + le £^ = 0. Although one can write an expression for 
rn(M, q, Q), it turns out to be more convenient to express q in terms of th, Q, M: 



r H / 8M 3 



2VM V ( 2 ^ 2 + Mr H 



2 \2 



1 . (3.11) 



As usual, the constants $ , <I> a in the expressions of A t (r),f? t (r) are found by imposing the regularity 
of the one-forms A, B on the horizon, which implies 



MQ j>° = Mq (3 12) 

2Q 2 + Mr 2 H \q- 1 Mr 2 ,' 



Thus, physically they correspond to the two chemical potentials associated with the system. 

A straightforward computation leads to the following expressions for the mass A4, electric charges Qa 
and Q a for the SU{2) and C7(l) fields, entropy S and Hawking temperature T#: 

M = ^-V, Q A = ^QV, Qa = L qV , S = \V2Mr H V, 
l07r 2tt 2tt 2 

1 8Q 6 + 8M 4 r 2 H + l2MQ 4 r 2 H + 6M 2 Q 2 r 4 H + M 3 (-16Q 2 + r 6 H ) 
H ~ 32V2n MV 2 r H (2Q 2 + Mr 2 H ) ' 

with V = J d 3 x; however, for the rest of this work, to simplify the expressions, we set V = 1, i.e. we shall 
work with mass, entropy and electric charge densities. 

A straightforward computation shows that the solutions satisfy the first law of thermodynamics, dAA = 
T H dS + $ A dQ A + $> a dQ a , and the Smarr law, M = 1{T H S + $ A Q A + <P a Q a ). 

3.2 Thermodynamic properties 

These U(l) x U(l) solutions possess a relatively complicated thermodynamics. Restricting for simplicity to 
a canonical ensemble, we study black branes holding the temperature Th, and the charges Qa, Qa fixed. 
The associated thermodynamic potential is the Helmholz free energy F[Tjj; Qa, Qa] = M — ThS . Thermo- 
dynamic stability requires the positivity of the specific heat at constant electric charges, C = Tff(dS/dTff). 

A useful relation here is T H = s ^/s ^+32)2 ^2+32^2/3 , which defines S(T H ; Q A , Q a ), although an 

explicit formula cannot be written in the general case. 

Analytic results are found only when discussing the limiting cases with a vanishing Q or q. The properties 
of the solutions with Q — (i.e. a consistent truncation of the model with a Z7 (1) field only, F$ = 0) arc 
discussed at length in [3] . No extremal configurations are found in this case, since the temperature is bounded 
from below, j 1 ^"™) > q p or am7 given Th > Tjy , there are two branches of solutions, one of them being 
thermally stable. 

By constrast, the solutions with q = (SU(2) gauge fields only, J M „ = 0), admit an extremal limit which 
is approached for M = Q 4 / 3 /2 1 / 3 . The entropy of the extremal solutions vanishes, a number of invariant 
quantities diverging in that limit. For nonextremal configurations one finds a single branch of solutions, with 



S = 4ttT h Q 



2/3 



32tt 2 T= , / ^ 32768tt 6 T^ _ L 65536^?^ \ / ^ 32768tt 6 T^ L 65536ti- 6 T6 



3Q 2 / 3 V 27Q 2 V 2 7Q 2 / V 27< 2 2 V 27 Q 



1/3 / , \ 1/3 



which possesses a positive specific heat. 

The solutions with two U(l) charges exhibit a complicated picture, which is governed by the value of the 
relative ratio q/Q. The picture in Figure 1 appears to be generic: for any Q ^ 0, the solutions with small 
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Figure 1: The entropy and free energy of the Abelian black brane solutions is shown for several values of the electric 
charges. 



enough q are thermally stable, the entropy increasing with the temperature (note that S(Th = 0) =/= 0, the 
geometry remaining regular in this limit). 

However, when increasing q we notice the occurance of three branches of solutions for some intermediate 
range of Tg. The physically relevant branch (which has less free energy) is the third one, which continues to 
T H — > oo (the large temperature behaviour is S = 1287r 3 Tj^ + 0{T H )). Also, the second branch is unstable 
since it possess a negative specific heat. 

For a more systematic discussion of the properties of the generic Abelian solutions, it turns out convenient 
to work with the following scaled quantities^] 

QA Qa S T H A $ A a <S> a 

qA = Cl M*7*' qa = Ci MV^ S = C2 M^> tH = C3 Mm> v = C5 M^> v = C6 M^> (3 ' 13) 

with Cl = ^1(| )V*, C2 = 5^(1)3/4 c 3 = 4 x 6 1 /4 7r 3/4 ! C4 = Cg = (4)1/* and c 6 = 2^(1)1/4 

being constant factors which have been chosen such that the expressions below take a simpler form. 
All the relevant quantities can then be expressed in terms of qA and s only: 



It is clear that all solutions satisfy the condition 

s 2 < i - <&, 

the upper bound being approached for solutions with SU(2) fields only. Moreover, the condition tg > 
imposses a lower bound for the reduced entropy: 



s 2 > U(q A ) - q 2 A , (3.15) 



where 



V3 / n \ 1/3 



2 The U(l) X {7(1) exact solution has an extra scaling symmetry M — > A 4 .M, Q fe — > A 3 Q fc , T H — > XT H , S — > A 3 S and 
<J> fc — > A$ fc , with k = (a, A) and A > an arbitrary constant. The quantities in (13. 1 3 P > arc left invariant by this transformation. 
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One can also show that the scaled £7(1) charge q a has a finite range, with 

Solutions with a maximal value of q a correspond to extremal black holes, with tjj = and s 2 = U(qA) — q\- 
From p,14p . the reduced entropy of the extremal solutions can also be written as 

1 

4 



g2 = l \-Ql + QaJql + Sq 2 A ) , (3.18) 



which is a nonvanishing quantity for q& ^ 0. 



4 The superfluid phase 

It is clear that the £7(1) x £7(1) solutions should possess non-Abelian generalizations. These configurations 
are found when enlarging the SU(2) ansatz to include a nonzero magnetic potential such that the gauge 
potential = A t (r)S I3 dt is approached only asymptotically. 

Following previous work [5], [5] on pure Einstein- Yang-Mills (EYM) solutions with vector hair, we choose 
an SU(2) gauge fields ansatz possessing both electric and magnetic potentials, while the £7(1) ansatz is still 
purely electric: 

A {1) = w(r)S n dx + A t (r)S I3 dt, B = B t (r)dt. (4.19) 

Also, as before, the dilaton field will depend on the r-coordinate only. This leads to a diagonal energy- 
momentum tensor and thus it is consistent to choose a diagonal metric ansatz. 



4.1 Zero modes for the £7(1) x £7(1) black brane 

Before discussing the general solutions, it is instructive to consider the perturbative limit of the problem. 
Then w(r) is treated as a small perturbation around the £7(1) x £7(1) solutions, w(r) = eW(r). After 
substituting into the linearized YM equations, one finds that W(r) solves 

,. 1 K' f\ . A 2 H 2 K 
W" + — + J —)W' + -S W = 0. 4.20 

r K f f 

For Q = one finds the following exact solution of the above equation 

W(r) =c + Cl (log(l - Cff) - log(l + , (4.21) 

(where Cq,Ci are arbitrary constants). As one can see, this solution posseses an essential logarithmic singu- 
larity at the horizon and thus cannot be treated as a perturbation. Thus we conclude that only solutions 
which are charged with respect the SU{2) fields may posses an instability. 

Although for Q ^ the equation (|4.20p does not appear to be solvable in terms of known functions, one 
can construct approximate solutions near the horizon and at infinity. As r — > rji , the function W(r) behaves 
as W(r) = b + 0(i — ru ) 2 , while, for large r, the approximate form of W(r) is W(r) = J/r 2 + 0(l/r 4 ), with 
b and J free parameters. Solutions interpolating between these asymptotics are constructed numerically 

The mechanism triggering the instability is similar to the pure EYM- A case [2], the magnetic gauge 
potential acquiring a tachyonic mass term near the horizon. Interestingly, the picture found for q = is 
rather similar to that valid for Q = since in this case too no solutions of (|4.20l) with correct asymptotics 
are found. We conclude that, somehow unexpectedly, both electric charges (associated with the SU(2) and 
£7(1) fields) should be nonvanishing for the existence of a normalizable zero mode. 



3 In this work we restrict our study to solutions with a monotonic behaviour of W(r). 
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Figure 2: Critical curves in the parameter space where static linear normalizable non-Abelian perturbations arise. 



Some results of the numerical integration are shown in Figure 2. There, the part of the parameter space 
above the curves corresponds to the unbroken phase, where only Abelian solutions exist. In Figure 2 (left) 
we show the critical curve in the {qA,q a ) plane corresponding to configurations unstable with respect to 
non-Abelian perturbations. One can see that the reduced SU(2) charge qA = ciQa/A4 3 ^ 4 has a finite range, 
< qA < 0.618 , an extremal configuration (with Tjf* — > 0) being approached for the maximal value of qA 
and q a — > 0.629 (corresponding to <fr a /<fr A ~ 0.704). In Figure 2 (right) we show the critical temperature 
as a function of the U(l) chemical potential $ a (both quantities are normalized w.r.t. the SU(2) chemical 
potential $ j4 ); note that the critical temperature is monotonically decreasing as we increase the ratio & a /$ A . 



4.2 Black holes with non-Abelian hair 
4.2.1 The equations and global charges 

The instability of the U(l) x £7(1) solution pointed out in the previous section can be viewed as an indication 
of the existence of a branch of non-Abelian solutions with nontrivial magnetic non-Abelian fields outside the 
horizon. 

In the numerical construction of these solutions, we adopt the following metric ansat^j], which was first 
proposed in [8] for the case of the pure EYM-A system: 

^ = W) + ^ iMr) + f2{r){dy2 + dz2) ) N ( r ) a2 ( r ) dt2 i with N(r) = + (4.22) 

Plugging (|4.22l) and (14.19)) into (j2~3"l) - (|2~4"f results in the equations of motiorQ 



3r 3 Nf 2 re 2 ^ f , A n r 2 A' 2 f 4 A 2 w 2 
2/ 2 9ym \ a Nc 



-4a<p 3 R /2 i 2 

+ 2a 2 + r 3N( t>' 2 + -^r 3 (3g 2 YM - e~ 2 ^ - 2e a % (4.23) 
o> _ 2r (if 2 a 1 e 2a * fi a , A 2 w 2 f 4 



° 3 V P T ' 9 2 ym r 2 w N 2 a 2 



4 The line element 1 13. 8 I t of the £7(1) X !7(1) solution can also be written in the form 114.221 1 (with /(r) = 1) by defining a new 
radial coordinate. However, this results in rather complicated expressions. 

5 One can see that the Chern-Simons term in 112.11 1 does not contribute to the equations of motion so that the gauge fields 
do not interact directly. 
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f 2p 2ar t> A 2 w 2 f' 2 

(r 3 Na<j)'y = aro 



9ym 



J2AI2 fi A2 2\ o„-4a0 2 
f^Nni 12 - ' - i * I 4- r 2 R' 2 + SyM fr" 2 '"*' - P a 4>V 2 

/ ivw ^ iVa 2 J + a 2 * + 8 1 ' 



A> „2a<p f4 r>/ 2a<p fi 

(e 2a0 r 3fty = ( e - 4a V^)' = 0, (e 2a *Nf A raw')' = - e —fl-A 2 w . 

o No o No 

These equations possess the following scaling symmetries (invariant functions are not shown) 

w 

(i)o^\o, A t — !> XA t , B t ^\B u (ii) f ->■ A/, w -)• (4.24) 

(m) r — > Ar, <?kjw -4 ^ym^ ^ _^ t ^ ^ _^ ^ ^ _^ A 4 rn, A 4 — > AA t , _B f — > A£?t, to — > Aw, 

A A 

with A > an arbitrary number. The symmetries (i) and (ii) are used to set o(oo) = 1, /(oo) = 1, while 
(in) is used to set gyM = 1 without any loss of generality, which fixes the AdS length scale, £ = 2y2- Note 
also that the last equation in (|4.23[) implies the existence of the first integral 

s^V^' (425) 

with q a constant fixing the U(l) electric charge. 

We consider again black branes with an horizon at r — r^, where N(r^) = 0. The non-extremal solutions 
have the following expansion as r — > r#: 

q 4 

= 4^ + °( r - r ff)' <r(r)=<T h + 0(r-r H ), f(r) = f h + 0(r-r H ) 2 , cj>(r) = O + 0(r - r H ), 

w(r) = w h + 0(r - r H ) 2 , A t (r) =V x (r- r H ) + 0(r - r H ) 2 , B t (r) = u x (r - r H ) + 0(r - r H ) 2 , (4.26) 

with the independent parameters {oh, fh, $0, w h, v\, Vi} which fix the coefficients of all higher order 
terms. 

We are interested in solutions approaching the U(l) x U(l) configurations asymptotically. We assume^ 
that, as r — > oo, w(r) vanishes and <j)(r) decays as 1/r 2 . A systematic analysis then reveals the following 
expansion of the solutions at large r: 



m(r) = M + 0(l/r 2 ) 1 o(r) = 1 - -— + 0(l/r 4 ), <f>(r) = - + 0(l/r 4 ), f(r) = 1 + J ± + 0(l/r% 
w (r) = ^ + 0(l/r 4 ), A t (r)=Z> A -%, B t (r) = $ a -\, (4.27) 



with {M, J, Q, q, D, f 4 , $ a } arbitrary coefficients. 

All physical quantities are fixed by the data at the horizon and at infinity. As in the Abelian case, the 
global charges are the mass and the SU(2) and U(l) electric charges, witfQ 

M = ±(M+%.) t Qa = ±Q, Q a = ± q , (4.28) 

while & A , $° are chemical potentials associated with the two gauge fields. The entropy and Hawking 
temperature of the solutions are given by 



c 1 3 rj, °h 



-(-(2e a *° + e - 2a *°) - 2 * ) 



rtJ 1 Op2cup 2 S p 4a0o„2T 

'H , J- , n _ad> n , _-2a6n\ Ze u l \ oe 1 



3 ^ ' o 2 h ' 3rl 



(4.29) 



3 The generic solutions have a more complicated asymptotic behaviour, with ui(r) = ujq — l 2 ti?^} 1 °^ r + . . . , 4>(r) 



fi 1 °^ r + . . . , which implies the existence of log terms also in the expression of the metric functions, e.g. m(r) = M + j^/3(/3 — 

4a) logr - ^- log 2 r + logr + ... . 

7 Note that, different from the pure EYM— A case, the total mass is not given by the asymptotic value of m(r), aquiring a 
contribution from the scalar field. 
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Figure 3: The profile of a typical non-Abelian solution is shown as a function of the radial coordinate. 

For completness, we mention that the boundary stress-tensor Tj as defined by (|2.6[) is diagonal, with the 
nonzero components: 

T * = ^37 (i?2 ~ 12/4 + M£2) ' jy y = JZ * = -kU D2 + Qh + Mei Jtt = ~^I iD2+Ml2) ' (4 ' 30) 

such that T\ = 0. 

4.2.2 Numerical solutions 

Eqs. (|4.23j) with boundary conditions (|4.26|) and (|4.27|) . respectively, have been solved numerically using 
a standard shooting method. In addition to using this algorithm, some solutions were also constructed by 
employing a collocation method for boundary- value ordinary differential equations equipped with an adaptive 
mesh selection procedure. We have confirmed that there is good agreement between the results obtained 
with these two different methods. 

As expected, some basic properties of these black branes are rather similar to those found in [8], [9] in 
the case of the purely EYM-A model. However, the solutions in the present work feature a second control 
parameter, which is the U(V) electric charge q a (or equivalently, the chemical potential $ a ). 

For all solutions, the functions a(r) and A t (r), B t (r) always increase monotonically with growing r. 
However, m(r), f{r), 4>(r) and w(r) may feature a more complicated behaviour, with local extrema. For 
sufficiently small u>h, all field variables remain close to their values for the Abelian configuration with the 
same (rn,Q,q)- Significant differences occur for large enough values of Wh and the effect of the magnetic 
fields on the geometry becomes increasingly more pronounced. The profiles of a typical solution illustrating 
these features are presented in Figure 3. 

In the numerical approach, we make use of the existence of the first integral (|4.25[) to fix the value of 
the electric charge associated with the U{\) field, which implies v\ = 8e 4a ^ h qah/rjj in the near horizon 
expansion (|4.26p . The scaling symmetry (iv) in (|4.24|) is used to set rg = 1, such that the only remaining 
control parameters are w(ru) and V±. 

We have studied in a systematic way families of solutions with fixed values of q between 0.5 and 7, the 
following picture being generic. First, the behaviour of solutions for arbitrary data on the horizon is such 
that at large r one finds w —> wo ^ and <f>{r) — > log r/r 2 (in which case the total mass as defined according 
to (|2.7|) diverges), or else there is a singularity at finite r. Given (u^, g; ry), solutions with the correct 
asymptotic behavioujfl are found only for a discrete set of values of (Vi, (f>o). Also, all solutions possess a 

8 These solutions here are also indexed by the node number of the magnetic potential w(r). It turns out that the configurations 
with nodes represent excited states whose energy is always greater than the energy of the corresponding nodeless configurations, 
and are therefore ignored in what follows. 
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Figure 4: Various parameters of the non-Abelian solutions are shown for several values of of the scale-free ratio 

7 = Qa/S. 

non-vanishing electric Qa associated with the SU(2) field. 

Moreover, for fixed (rn,q), one finds a branch of non-Abelian solutions for < Wh < w h nax . Along this 
branch, the Hawking temperature decreases, an extremal configuration being approached^ for the maximal 
value of Wh ■ The numerical construction of the solutions with Tjj = requires a different metric ansatz than 
(|4.22p and is beyond the scope of this work. However, based on the results in the near-extremal case, we 
expect the extremal solutions to share the basic properties of their general Abelian counterparts, possesing 
a regular horizon with non- vanishing entropy. 

Some results of the numerical integration are shown in Figure 4. There we employ scale-free quantities 
defined in ()3.13j) . which are invariant under the scaling transformation (iv) in (|4.24j) : also, we have found 
it convenient to define 7 = Q a /S as a second scale-free control parameter. One can see that, for all 
cases we considered, non-Abelian solutions exist only for values of the Hawking temperature smaller than 
a critical temperature This is the temperature at which the £7(1) x £7(1) solution admits a static 

linearised perturbation, with nonvanishing but infinitesimally small w. Moreover, the dependece of the 
order parameter J on the Hawking temperature is similar to that found in the literature for the 7 = case 
(i.e. a EYM-A model). Also, as expected, we have found that the difference in the free energy density, F, 
between a non-Abelian solution and the £7(1) x £7(1) solution with the same temperature and electric charges 
is negative, and thus the non-Abelian solution is thermodynamically favoured. 



5 Further remarks. 

In this paper we have studied electrically charged black branes of the TV = 4+ SU(2) x £7(1) gauged 
supergravity model with AdSs asymptotics. Apart from the Abelian £7(1) x £7(1) configurations, we have 
given numerical evidence that this model possesses also solutions with a non-vanishing magnetic SU(2) 
fields. Remarkably, these emerge as perturbations of the Abelian configurations at some finite temperature 
depending on the values of the electric charges, and can be viewed as p— wave superfluids. Moreover, by 
using the relations in [TU] , [U , one can uplift these configurations to ten dimensional type IIB supergravity 
and D = 11 supergravity. This provides an explicit stringy construction of holographic superfluids. 

Our study should be viewed only as a preliminary investigation of the simplest non-Abelian solutions of 
the Af = 4 + model featuring superfluid properties. Various properties of these black branes remain to be 

9 Our numerical code usually provided good quality solutions for Th ^ T^'/10. 
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investigated. For example, it would be interesting to compute the conductivity of the solutions or to explore 
the connection with the unbalanced mixtures discussed recently in |15j . 

Moreover, we expect the AT — 4 + model to possess a variety of other electrically charged black brane 
solution^. They would be found for a different (typically more complicated) matter field ansatz than (14.191) . 
In this context, it is interesting to study solutions in which the Chern-Simons term enters the dynamics n l. 
To this end, we have considered non-Abelian black branes possesing a purely magnetic SU (2) field, with 

A {1) = w(r)(6 n dx + S I2 dy + S I3 dz), (5.31) 

and an electric U(l) field, B = B t (r)dt. This leads to an isotropic energy-momentum tensor, T* — Ty — T*, 
in which case a suitable metric Ansatz is given by (|4.22[) with f(r) — 1. Then (|2.3[) - (|2.4j) yield five equations 
of motion for m, a, w, B t and <j> which were solved numerically. Our results show that the properties of 
these solutions differ substantially from those found in the anisotropic case, discussed in Section 4. First, 
when treating w(r) as a perturbation around the Abelian solution (which is (|3.8[) - (I3.10I) with Q = 0), the 
linearized equation can be solved in closed form. However, the solution looks very similar to (|4.2ip (with a 
log(r — rjj) term), with the result that no normalizable zero mode is found. Also, different from the case 
of anisotropic non-Abelian black branes, we could not find non-perturbative solutions with w(r) — > as 
r — > oo. As a result, the mass of the solutions computed according to the counterterm prescription given in 
Section 2, diver gea 12 l. We therefore conclude that these isotropic black brane non-Abelian solutions cannot 
be interpreted as holographic superfluids. 

However, the situation is likely to be different for a more general case featuring an anisotropic SU(2) 
field and a purely magnetic U(l) field (thus beyond the simple Ansatz (14.19p ~). Superconducting black brane 
solutions of this type have been studied recently in a truncation of Af = 4 + Romans' model with a vanishing 
dilaton and an arbitrary Chern-Simons coupling constant |16) . 

We hope to return elsewhere with a systematic study of these aspects. 
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